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Abstract We review the physics of relativistic shocks, which are often in¬ 
voked as the sources of non-thermal particles in pulsar wind nebulae (PWNe), 
gamma-ray bursts (GRBs), and active galactic nuclei (AGN) jets, and as possi¬ 
ble sources of ultra-high energy cosnric-rays. We focus on particle acceleration 
and magnetic field generation, and describe the recent progress in the field 
driven by theory advances and by the rapid development of particle-in-cell 
(PIC) simulations. In weakly magnetized or quasi parallel-shocks ( i.e. where 
the magnetic field is nearly aligned with the flow), particle acceleration is 
efficient. The accelerated particles stream ahead of the shock, where they gen¬ 
erate strong magnetic waves which in turn scatter the particles back and forth 
across the shock, mediating their acceleration. In contrast, in strongly magne¬ 
tized quasi-perpendicular shocks, the efficiencies of both particle acceleration 
and magnetic field generation are suppressed. Particle acceleration, when effi¬ 
cient, modifies the turbulence around the shock on a long time scale, and the 
accelerated particles have a characteristic energy spectral index of s 7 ~ 2.2 
in the ultra-relativistic limit. We discuss how this novel understanding of par¬ 
ticle acceleration and magnetic field generation in relativistic shocks can be 
applied to high-energy astrophysical phenomena, with an emphasis on PWNe 
and GRB afterglows. 
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1 Introduction 

In pulsar wind nebulae (PWNe), gamma-ray bursts (GRBs), and jets from ac¬ 
tive galactic nuclei (AGNs), signatures of non-thermal processes are revealed 
by power-law radiation spectra spanning an extremely wide range of wave¬ 
lengths, from radio to X-rays, and beyond. Yet, it is still a mystery how the 
emitting particles can be accelerated up to ultra-relativistic energies and how 
the strong magnetic fields are generated, as required in order to explain the 
observations. In most models, non-thermal particles and near-equipartition 
fields are thought to be produced at relativistic shock fronts, but the details 
of the mechanisms of particle acceleration and magnetic field generation are 
still not well understood. 

Particle acceleration in shocks is usually attributed to the Fermi process, 
where particles are energized by bouncing back and forth across the shock. 
Despite its importance, the Fermi process is still not understood from first 
principles. The highly nonlinear coupling between accelerated particles and 
magnetic turbulence - which is generated by the particles, and at the same 
time governs their acceleration - is extremely hard to incorporate in analytic 
models. Only in recent years, thanks to major breakthroughs on analytical and 
numerical grounds, has our understanding of the Fermi process in relativistic 
shocks significantly advanced. This is the subject of the present review. 

Relativistic shocks pose some unique challenges with respect to their non- 
relativistic counterparts. For example, the distribution of accelerated particles 
can no longer be approximated as isotropic if the shock is relativistic. In a 
relativistic shock, the electric and magnetic fields significantly mix as one 
switches between upstream and downstream frames of reference. And unlike 
non-relativistic shocks, where some aspects of the theory can be tested by 
direct spacecraft measurements, relativistic shocks are only constrained by 
remote observations. For recent reviews of relativistic shocks, see Bykov and 
Treumann (2011); Bykov et al. (2012). 

This chapter is organized as follows. First, we review recent analytical 
advances on the theory of particle acceleration in relativistic shocks, arguing 
that the accelerated particle spectrum and its power-law slope in the ultra- 
relativistic limit, s 7 = — dlog N/dlogj ~ 2.2 (where 7 is the particle Lorentz 
factor), are fairly robust (Section 2). Here, we assume a priori that some 
magnetic turbulence exists on both sides of the shock, such that the Fermi 
process can operate. Next, we describe the plasma instabilities that are most 
relevant for generating this turbulence (Section 3), stressing the parameter 
regime where the so-called Weibel (or “filamentation”) instability - which is 
often thought to mediate the Fermi process in weakly magnetized relativistic 
shocks - can grow. Then, we summarize recent findings from particle-in-cell 
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(PIC) simulations of relativistic shocks, where the non-linear coupling between 
particles and magnetic waves can be captured from first principles (Section 
4). Finally, we describe the astrophysical implications of these results for the 
acceleration of ultra high energy cosmic rays (UHECRs) and for the radiative 
signatures of PWNe and GRB afterglows (Section 5; for a review of PWNe, 
see Kargaltsev et al. (2015) in the present volume; for a review of GRBs, see 
Racusin et al. (2015) in the present volume). We briefly conclude in Section 6. 


2 Particle Acceleration in Relativistic Shocks 

Diffusive (Fermi) acceleration of charged particles in collisionless shocks is 
believed to be responsible for the production of non-thermal distributions of 
energetic particles in many astronomical systems (Blandford and Eichler 1987; 
Maikov and Drury 2001, but see, e.g. Arons and Tavani (1994) for a discussion 
of alternative shock acceleration processes). The Fermi acceleration process in 
shocks is still not understood from first principles: particle scattering in colli¬ 
sionless shocks is due to electromagnetic waves formed around the shock, but 
no present analytical formalism self-consistently calculates the generation of 
these waves, the scattering and acceleration of particles, and the backreaction 
of these particles on the waves and on the shock itself. 

The theory of particle acceleration was first developed mainly by evolv¬ 
ing the particle distribution under some Ansatz for the scattering mechanism 
(e.g. diffusion in pitch angle), within the “test particle” approximation, where 
modifications of wave and shock properties due to the high energy particles 
are neglected. This phenomenological approach proved successful in explaining 
the spectrum of relativistic particle distributions inferred from observations, 
although a more careful approach is needed to account for the energy fraction 
deposited in each particle species (electrons, positrons, protons, and possibly 
heavier ions), and to test the Ansatz of the scattering prescription. 

For non-relativistic shocks, the linear theory of diffusive particle acceler¬ 
ation, first developed in 1977 (Krymskii 1977; Axford et al. 1978; Bell 1978; 
Blandford and Ostriker 1978), yields a power-law distribution d 3 N/d 3 p oc p~ Sp 
of particle momenta p , with a spectral index 

s p = s 7 + 2 = 3/3 u /(j3 u — pd). (1) 

Here, /3 is the fluid velocity normalized to the speed of light c in the frame of 
the shock, which is assumed planar and infinite, and subscripts u (d) denote the 
upstream (downstream) plasma. For strong shocks in an ideal gas of adiabatic 
index r = 5/3, this implies s p = 4 (i.e. s 7 = 2; constant energy per logarithmic 
energy interval, since p 2 d 3 N/d 3 p oc p~ 2 ), in agreement with observations. 

The lack of a characteristic momentum scale, under the above assumptions, 
implies that the spectrum remains a power-law in the relativistic case, as 
verified numerically (Bednarz and Ostrowski 1998; Achterberg et al. 2001). 
The particle drift downstream of the shock implies that more particles are 
moving downstream than upstream; this anisotropy is of order of p u when 
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measured in the downstream frame (Keshet and Waxman 2005). Thus, while 
particle anisotropy is negligible for non-relativistic shocks, the distribution 
becomes highly anisotropic in the relativistic case, even when measured in 
the more isotropic downstream frame. Consequently, one must simultaneously 
determine the spectrum and the angular distribution of the particles, which is 
the main difficulty underlying the analysis of test particle acceleration when 
the shock is relativistic. 

Observations of GRB afterglows led to the conclusion that highly rela¬ 
tivistic collisionless shocks produce a power-law distribution of high energy 
particles with s p = 4.2 ± 0.2 (Waxman 1997; Freedman and Waxman 2001; 
Berger et al. 2003). This triggered a numerical investigation of particle accel¬ 
eration in such shocks, showing that s p indeed approaches the value of 4.2 for 
large shock Lorentz factors ( 7 „ = (1 — p^) 1 ^ 2 ^ 1); i n agreement with GRB 
observations, provided that particle scattering is sufficiently isotropic. 

The spectral index s p was calculated under the test particle approximation 
for a wide range of shock velocities, various equations of state, and different 
scattering prescriptions. This was achieved by approximately matching nu¬ 
merical eigenfunctions of the transport equation between upstream and down¬ 
stream (Kirk and Schneider 1987; Heavens and Drury 1988; Kirk et al. 2000), 
by Monte Carlo simulations (Bednarz and Ostrowski 1998; Achterberg et al. 
2001; Ellison and Double 2002; Lemoine and Pelletier 2003; Niemiec and Os¬ 
trowski 2004; Lemoine and Revenu 2006; Ellison et al. 2013), by expanding 
the distribution parallel to the shock front (Keshet and Waxman 2005), and 
by solving for moments of the angular distribution (Keshet 2006). 

These studies have assumed rest frame diffusion in pitch angle or in the 
angle between particle velocity and shock normal. These two assumptions yield 
similar spectra in the limit of ultra-relativistic shocks (Ostrowski and Bednarz 
2002 ). As discussed later in this review, one expects these assumptions to hold 
at relativistic shocks. However, some scenarios involve the conversion of the 
accelerated species into a neutral state and then back - e.g. proton to neutron 
and then back to proton via photo-hadronic interactions (Derishev et al. 2003) 
or electron to photon and then back to electron through Compton and pair 
production interactions (Stern and Poutanen 2008) in which case the particle 
may have time to suffer a large angle deflection upstream of the shock, leading 
to large energy gains and generically hard spectra (Bednarz and Ostrowski 
1998; Meli and Quenby 2003; Blasi and Vietri 2005). 

For isotropic, small-angle scattering in the fluid frame, expanding the par¬ 
ticle distribution about the shock grazing angle (Keshet and Waxman 2005) 
leads to a generalization of the non-relativistic Eq. (1) that reads 


s p — (3/3u — 2PuPd, + Pd)/(Pu — Pd) > 


( 2 ) 


in agreement with numerical studies (Kirk et al. 2000; Achterberg et al. 2001) 
over the entire range of p u and Pd- In particular, in the ultra-relativistic shock 
limit, the spectral index becomes 


s p {p u -> 1 , p d -> 1 / 3 ) = 38/9 = 4.222 ... 


( 3 ) 
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The spectrum is shown in Fig. 1 for different equations of state, as a function 
of the shock four-velocity 7„/3 u . 



Fig. 1 Spectral index according to Eq. (2) (Keshet and Waxman 2005, curves) and to 
a numerical eigenfunction method (Kirk et al. 2000, symbols), as a function of , for 
three different types of shocks (Kirk and Duffy 1999): a strong shock with the Jiittner/Synge 
equation of state (solid curve and crosses), a strong shock with fixed adiabatic index F = 4/3 
(dashed curve and x-marks), and for a relativistic gas where /3 u /b = 1/3 (dash-dotted curve 
and circles). 


The above analyses assumed that the waves scattering the particles move, 
on average, with the bulk fluid velocity. More accurately, one should replace 
/3 by the mean velocity of the waves that are scattering the particles. In the 
shock precursor (see §3.1), the scattering waves are expected to be slower than 
the incoming flow, leading to a softer spectrum (smaller f3 u in Eq. 2). 

Small-angle scattering can be parameterized by the angular diffusion func¬ 
tion V = {{AO) 2 /At), where 0 is the angle of the particle velocity, taken here 
with respect to the shock normal, and angular brackets denote an ensemble 
average. The function V = V{0,p,z) should be specified on both sides of the 
shock, and in general depends on 0 , on the particle momentum p, and on its 
distance 2 from the shock front. 

For scattering off waves with a small coherence length A -C rz,, where 
77 = ( pc/eB ) is the Larmor radius, roughly ( 77 ,/A) 2 uncorrelated scatter¬ 
ing events are needed in order to produce an appreciable deflection, so T> ~ 
r\c/ A oc p 2 (Achterberg et al. 2001; Pelletier et al. 2009). Here, B is the 
magnetic field, and e is the electron’s charge. Simulations (Sironi et al. 2013) 
confirm this scaling at early times; some implications are discussed in Section 
4. The precise dependence of V upon 2 is not well known. It is thought that 
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T> slowly and monotonically declines away from the shock, as the energy in 
self-generated fields decreases. However, the extents of the upstream precursor 
and downstream magnetized region are not well constrained observationally, 
and in general are numerically unaccessible in the foreseeable future. 

For an evolved magnetic configuration, it is natural to assume that the 
diffusion function is approximately separable in the form T> = D(9)D2(p, z). 
Here, D 2 (which may be approximately separable as well, but see Katz et al. 
2007) can be eliminated from the transport equation by rescaling z, such that 
the spectrum depends only on the angular part D{6). 


The spectrum is typically more sensitive to the downstream diffusion func¬ 
tion Dd than it is to the upstream D u . In general, an enhanced Dd along 
(opposite to) the flow yields a softer (harder) spectrum; the trend is roughly 
reversed for D u (Keshet 2006). Thus, the spectrum may deviate significantly 
from that of isotropic diffusion, in particular in the ultra-relativistic limit (Kirk 
et al. 2000; Keshet 2006). However, the spectral slope s is not sensitive to lo¬ 
calized changes in D at angles perpendicular to the flow (Keshet 2006). For 
roughly forward-backward symmetric scattering in the downstream frame, as 
suggested by PIC simulations, s is approximately given by its isotropic diffu¬ 
sion value in Eq. 2 (Keshet et al., in preparation). 

Particle acceleration is thought to be efficient, at least in weakly magnetized 
or quasi-parallel shocks, as discussed below. Thus, the relativistic particles are 
expected not only to generate waves, but also to slow down and heat the 
bulk plasma (Blandford and Eichler 1987). As particles with higher energies 
are expected to diffuse farther upstream and slow the plasma, lower-energy 
particles are effectively accelerated by a slower upstream. Consequently, if the 
scattering waves are assumed to move with the bulk plasma, the spectrum 
would no longer be a power-law. However, this effect may be significant only 
for mildly relativistic shocks, with Lorentz factors below 7 ^ ~ 3 (Ellison and 
Double 2002; Ellison et al. 2013). 

To understand the energy, composition, and additional features of the ac¬ 
celerated particles, such as the acceleration time and energy cutoffs, one must 
not only analyze the scattering of these particles (for example, by deriving V), 
but also address the injection problem, namely the process by which a small 
fraction of particles becomes subject to significant acceleration. Such effects 
were investigated using Monte-Carlo techniques (Bednarz and Ostrowski 1998; 
Ellison et al. 2013), in the so-called “thermal leakage” model, where fast parti¬ 
cles belonging to the downstream Maxwellian are assumed to be able to cross 
the shock into the upstream. More self-consistent results on particle injection 
based on PIC simulations are presented in Section 4. To uncover the physics 
behind the injection and acceleration processes, we next review the generation 
of electromagnetic waves in relativistic shocks. 
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3 Plasma Instabilities in Relativistic Shocks 

3.1 The Shock Precursor 

The collisionless shock transition is associated with the build-up of some elec¬ 
tromagnetic barrier, which is able to slow down and nearly isotropize the in¬ 
coming unshocked plasma. In media of substantial magnetization 1 , a > 10 -2 , 
this barrier can result from the compression of the background magnetic field 
(as a result of the Lorentz transformation to the frame of a relativistic shock, 
the most generic configuration is that of a quasi-perpendicular field), while 
at lower magnetizations, it is understood to arise from the generation of in¬ 
tense micro-turbulence in the shock “precursor”, as explained hereafter and 
illustrated in Fig. 2. 

At high magnetization, the gyration of the ambient particles in the back¬ 
ground compressed magnetic field can trigger a synchrotron maser instability, 
which sends precursor electromagnetic waves into the upstream (Langdon et al. 
1988; Hoshino and Arons 1991; Hoshino et al. 1992; Gallant et al. 1992). As 
incoming electrons and positrons interact with these waves, they undergo heat¬ 
ing (Hoshino 2008; Sironi and Spitkovsky 2011b), but acceleration seemingly 
remains inefficient (Section 4.1). 

At magnetizations a < 10~ 2 , the interpenetration of the incoming back¬ 
ground plasma and the supra-thermal particles, which have been reflected on 
the shock front or which are undergoing Fermi cycles around the shock, leads 
to anisotropic micro-instabilities over an extended region in front of the shock, 
called the “precursor” here. These instabilities then build up a magnetic bar¬ 
rier, up to a level 2 cb ~ 10 —2 —10 —1 , sufficient to deflect strongly the incoming 
particles and thus mediate the shock transition. This picture, first envisioned 
by Moiseev and Sagdeev (1963), has been recently demonstrated in ab initio 
PIC simulations (Spitkovsky 2005; Kato 2007; Spitkovsky 2008a). The gener¬ 
ation of micro-turbulence in the shock precursor is thus a key ingredient in 
the formation of the shock and in the development of the Fermi process, as 
anticipated analytically (Lemoine et al. 2006) and from Monte Carlo simula¬ 
tions (Niemiec et al. 2006), and demonstrated by PIC simulations (Spitkovsky 
2008b; Sironi and Spitkovsky 2009a, 2011b), see hereafter. 

As seen in the background plasma (upstream) rest frame, the supra-thermal 
particles form a highly focused beam, with an opening angle ~ l/ 7 „ and a 
mean Lorentz factor 7 |, u ~ y 2 . In contrast, boosting back to the shock frame, 
this supra-thermal particle distribution is now open over ~ 7r/2, with a mean 

1 The magnetization is defined as a = B 2 / [4 Try,,(7,, — l)n'mc 2 ] in terms of B, the 
large-scale background magnetic field measured in the shock front rest frame, and n' , the 
proper upstream particle density. The mass m is rn p for an electron-proton shock, and m e 
for an electron-positron shock, i.e. it corresponds to the mass of the particles which carry 
the bulk of the inertia. For a perpendicular shock, in which the background magnetic field 
in the shock frame is perpendicular to the flow, the magnetization can also be written as 
a = (u\ / c) 2 , with u\ the Alfven four-velocity of the upstream plasma. 

2 The parameter eg denotes the magnetization of the turbulence, eg = 
SB 2 / [47T7 u (7 u — l)n'mc 2 ], where SB is the fluctuating magnetic field. 
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Fig. 2 Phase diagram of relativistic collisionless shocks in the plane ( 7 u,cr); this figure 
assumes > 10 and £ C r = 0 . 1 , where the parameter £ C r = e cr / [ 7 u( 7 u — l)rc/rac 2 ] char¬ 
acterizes the energy density of supra-thermal particles (e cr ) relative to the incoming energy 
flux, as measured in the shock rest frame. In region 1, the shock transition is initiated by 
magnetic reflection in the compressed background field, while in regions 2 — 5, the magnetic 
barrier is associated to the growth of micro-instabilities, as indicated. The solid diagonal line 
indicates values of a and 7 u above which the filamentation instability would not have time 
to grow, in the absence of deceleration resulting from the compensation of the perpendicular 
current of the supra-thermal particles gyrating in the background field. See Section 3.1 and 
Lemoine et al. (2014a) for a detailed discussion. 


Lorentz factor 7 | sh > 7 ^, while the incoming plasma is highly focused, with 
a mean Lorentz factor j u . A host of micro-instabilities can in principle de¬ 
velop in such anisotropic configurations, see the general discussion by Bret 
(2009). However, in the deep relativistic regime, the restricted length scale of 
the precursor imposes a strong selection of potential instabilities, since a back¬ 
ground plasma volume element remains subject to the growth of instabilities 
only while it crosses the precursor. In the shock rest frame, this time scale is 
t x b — w ” 1 in the presence of a quasi-perpendicular background field 3 (a com¬ 
mon field geometry in relativistic flows), or t x , 5 B — 1 (w p A, 5 ,b/c) 1 

if the scattering is dominated by short scale turbulence of magnetization es 
and coherence length A sb (assuming that the waves are purely magnetic in 
the rest frame of the background plasma), see e.g. Milosavljevic and Nakar 


3 oj c = eB\ u /mc represents the upstream frame cyclotron frequency (and Bi u is the 

magnetic field in the upstream frame) while u) p = (A-nri!e~ /rn .) 1 1 denotes the plasma 
frequency. 
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(2006a), Lemoine and Pelletier (2010) and Plotnikov et al. (2013). This small 
length scale implies that only the fastest modes can grow, which limits the 
discussion to a few salient instabilities. 

Before proceeding further, one should stress that the above estimates for 
t x do not account for the influence of particles accelerated to higher energies, 
which can propagate farther into the upstream plasma and thus seed instabili¬ 
ties with smaller growth rate and on larger spatial scales. While such particles 
do not carry the bulk of the energy if the spectral index s 7 > 2, it is antici¬ 
pated that they should nevertheless influence the structure of the precursor, 
see in particular Milosavljevic and Nakar (2006a), Katz et al. (2007), Medvedev 
and Zakutnyaya (2009), Pelletier et al. (2009) and Reville and Bell (2014) for 
general analytical discussions, as well as Keshet et al. (2009) for an explicit 
numerical demonstration of their potential influence. Similarly, the above es¬ 
timates do not make a distinction between electron-positron and electron-ion 
shocks; in particular, it is understood that lo c and w p refer to the species which 
carries the bulk of the energy ( i.e. ions for electron-ion shocks). PIC simula¬ 
tions have demonstrated that in electron-ion shocks, electrons are heated in 
the precursor to nearly equipartition with the ions, meaning that in the shock 
transition their relativistic inertia becomes comparable to that of ions ( e.g. 
Sironi et al. 2013); hence one does not expect a strong difference between the 
physics of electron-positron and electron-ion shocks from the point of view of 
micro-instabilities, and unless otherwise noted, this difference will be omitted 
in the following. The microphysics of electron heating in the precursor never¬ 
theless remains an important open question, see Gedalin et al. (2008), Gedalin 
et al. (2012), Plotnikov et al. (2013) and Kumar et al. (2015) for recent discus¬ 
sions of this issue; indeed, the average Lorentz factor of electrons at the shock 
transition directly impacts the peak frequency of the synchrotron radiation of 
relativistic blast waves. 

In the context of relativistic weakly magnetized shocks, the most celebrated 
instability is the Weibel-like filamentation mode, which develops through a 
charge separation in the background plasma, triggered by magnetic fluctua¬ 
tions which segregate particles of opposite charges into current filaments of 
alternate polarity (e.g. Gruzinov and Waxman 1999; Medvedev and Loeb 
1999; Wiersma and Achterberg 2004; Lyubarsky and Eichler 2006; Achterberg 
and Wiersma 2007; Achterberg et al. 2007; Lemoine and Pelletier 2010; Bret 
et al. 2010; Rabinak et al. 2011; Lemoine and Pelletier 2011; Nakar et al. 2011; 
Shaisultanov et al. 2012). The current carried by the particles then positively 
feeds the magnetic fluctuations, leading to fast growth, even in the absence of a 
net large-scale magnetic field. In the rest frame of the background plasma, this 
instability grows in the linear regime as fast as 4 ~ £c/ 2 w p , with maximum 
growth on scales of the order of c/w p ; in the filament rest frame, this instability 
is mostly of magnetic nature, i.e. -Kw ~ 0. Several branches of this instability 
have been discussed in the literature, in particular the “oblique mode”, which 

4 The maximal growth rate of the Weibel instability is related to the plasma frequency of 

1/2 

the beam of supra-thermal particles, a; p b, though xsuj ~ cj p b, with o; p b — £ C r ^p- 
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involves a resonance with electrostatic modes. Even though this latter mode 
grows slightly faster than the fully transverse filamentation mode, it suffers 
from Landau damping once the electrons are heated to relativistic tempera¬ 
tures, while the transverse filamentation mode appears relatively insensitive 
to temperature effects. Thus, at a first order approximation, the transverse 
filamentation mode indeed appears to dominate the precursor at very low 
magnetizations. Its non-linear evolution, however, remains an open question; 
analytical estimates suggest that it should saturate at values eg <C 10 -2 via 
trapping of the particles (Wiersma and Achterberg 2004; Lyubarsky and Eich- 
ler 2006; Achterberg and Wiersma 2007; Achterberg et al. 2007), while PIC 
simulations see a continuous growth of magnetic energy density even when 
the non-linear filamentary structures have been formed ( e.g. Keshet et al. 
2009; Sironi et al. 2013). Whether additional instabilities such as a kinking of 
the filaments contribute in the non-linear phase thus remains debated, see for 
instance Milosavljevic and Nakar (2006b). 

At moderate magnetization levels, another fast instability can be triggered 
by the perpendicular current (transverse to both the magnetic field and the 
shock normal) seeded in the precursor by the supra-thermal particles during 
their gyration around the background field (Lemoine et al. 2014a, b). The com¬ 
pensation of this current by the background plasma on its entry into the pre¬ 
cursor leads to a deceleration of the flow, which modifies somewhat the effective 
timescale available for the growth of plasma instabilities, and destabilizes the 
modes of the background plasma. The growth rate for this instability can be 
as large as ~ w p , indicating that it can compete with the Weibel filamen¬ 
tation mode at moderate magnetizations. If the supra-thermal particle beam 
carries a net charge (in the shock rest frame), or a net transverse current, other 
similar instabilities are to be expected {e.g. Pelletier et al. 2009; Casse et al. 
2013; Reville and Bell 2014). The phase space study of Lemoine et al. (2014a) 
concludes that the filamentation mode likely dominates at magnetization lev- 
si 8 a 10 _ ', while this perpendicular current-driven instability dominates 
at 10~ 3 < cr < 10 —2 ; in between, both instabilities combine to form a com¬ 
plex precursor structure. Interestingly, these results do not seem to depend on 
the shock Lorentz factor, in good agreement with PIC simulations (Sironi and 
Spitkovsky 2009a, 2011b; Sironi et al. 2013). 

Finally, one should mention the particular case of quasi-parallel (sublumi¬ 
nal) configurations: there, a fraction of the particles can in principle escape 
to infinity along the magnetic field and seed other, larger scale, instabilities. 
One prime candidate is the relativistic generalization of the Bell streaming 
instability {e.g. Milosavljevic and Nakar 2006a; Reville et al. 2006), which 
is triggered by a net longitudinal current of supra-thermal particles; this in¬ 
stability has indeed been observed in PIC simulations (Sironi and Spitkovsky 
2011b). Of course, such a parallel configuration remains a special case in the 
deep relativistic regime. In mildly relativistic shock waves, with 7 U /3 U ~ 1, 
locally parallel configurations become more frequent, hence one could expect 
such instabilities to play a key role in seeding large scale turbulence. 
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3.2 Downstream Magnetized Turbulence 

How the magnetized turbulence evolves downstream of the shock is an impor¬ 
tant question, with direct connections to observations. The previous discussion 
suggests that the coherence length of the fields generated in Weibel-like insta¬ 
bilities should be comparable to the plasma skin-depth, c/ui p . However, mag¬ 
netic power on such scales is expected to decay rapidly through collisionless 
phase mixing (Gruzinov 2001), while modeling of GRB afterglow observations 
rather indicates that magnetic fields persist over scales ~ 10' —10 9 c/w p down¬ 
stream (Gruzinov and Waxman 1999). 

In a relativistic plasma, small-scale turbulence is dissipated at a damping 
rate 5 ~ —k 3 c 3 /ojp (Chang et al. 2008; Lemoine 2015) as a function of the 
wavenumber k, indicating that small scales are erased early on. Larger modes 
can survive longer; power on scales exceeding the Larmor radius of the bulk 
plasma decays on long, oc k 2 MHD scales (Keshet et al. 2009). It is not 
clear at present whether the small-scale turbulence manages to evolve to larger 
scales through inverse cascade effects ( e.g. Medvedev et al. 2005; Katz et al. 
2007; Zrake 2014), whether it is dissipated but at a rate which allows to match 
the observations (Lemoine 2013; Lemoine et al. 2013), or whether a large-scale 
field is seeded in the downstream plasma by some external instabilities (e.g. 
Sironi and Goodman 2007; Couch et al. 2008; Levinson 2009). 

PIC simulations have quantified the generation of upstream current fila¬ 
ments by pinching instabilities (e.g. Silva et al. 2003; Frederiksen et al. 2004; 
Jaroschek et al. 2005; Spitkovsky 2005, 2008a; Chang et al. 2008), and resolved 
the formation of shocks in two- and three-dimensional (2D and 3D) pair plasma 
(Spitkovsky 2005; Kato 2007; Chang et al. 2008; Keshet et al. 2009; Sironi and 
Spitkovsky 2009a; Haugbplle 2011; Sironi and Spitkovsky 2011b; Sironi et al. 
2013) and ion-electron plasma (Spitkovsky 2008a; Martins et al. 2009; Sironi 
et al. 2013). These simulations revealed a rapid decay of the magnetic field 
downstream at early times (Gruzinov 2001; Chang et al. 2008). Yet, a slow 
evolution of the plasma configuration takes place on > 10 3 /u;p timescales, 
involving a gradual increase in the scale of the magnetic structures, and con¬ 
sequently their slower dissipation downstream (Keshet et al. 2009). 

This long-term evolution is driven entirely by the high-energy particles ac¬ 
celerated in the shock; it is seen both upstream (e.g. in the precursor) and 
downstream, both of which become magnetized at increasingly large distances 
from the shock, and with an increasingly flat magnetic power-spectrum down¬ 
stream (Keshet et al. 2009). A flatter magnetic power spectrum at the shock 
implies a larger fraction of the magnetic energy stored in long-wavelength 
modes, which may survive farther from the shock. Indeed, the index of a power- 
law spectrum of magnetic fluctuations directly controls how fast the magnetic 
energy density, integrated over wavenumbers, decays behind the shock (Chang 


5 The shock crossing conditions imply that the relativistic plasma frequency of the shocked 
downstream plasma is roughly the same as the plasma frequency of the upstream plasma; 
no distinction will be made here between these quantities. 
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Fig. 3 Pair plasma evolution within 1000 c/lj p of the shock. The simulation is performed in 
the downstream frame, and the upstream flow moves with a Lorentz factor 7 r = 15 (so, 7 r is 
the relative Lorentz factor between the upstream and downstream regions). The normalized 
transverse magnetic field sign(13) eg (color scale stretched in proportion to 4 to highlight 
weak features) is shown at (a) early (ii = 2250 cjp 1 ), and (b) late (£2 = 11925 cj^ 1 ) times. 
Here Ax = x — x s h is the distance from the shock, with x s h (dashed vertical line) defined 
as the location of median density between far upstream and far downstream. Also shown 
are the transverse averages (at t±, dashed blue, and £ 2 ? solid red) of (c) the electromagnetic 
energy eem = [(H 2 + _E 2 )/ 87 r]/[( 7 r — l)^y r n'mc 2 ] (with E the electric field amplitude in the 
downstream frame, included in the definition of €em because in the simulation frame the 
induced electric field in the upstream medium is E ~ B) normalized to the upstream kinetic 
energy, (d) density normalized to the far upstream density n u = 7 ' r n f , and (e) particle 
momentum 7 /? (with (3 the velocity in c units) in the x-direction averaged over all particles 
(higher ( 7 (3 X )) and over downstream-headed particles only. 


et al. 2008; Lemoine 2015); the scale-free limit corresponds to a flat magnetic 
power spectrum (Katz et al. 2007). 

Properly capturing the backreaction of high energy particles requires large 
simulation boxes and large particle numbers, to guarantee that the largest scale 
fields and the highest energy particles are included. The largest available sim¬ 
ulations at the present, with length L and time T scales of (Tojp/c) 2 (Tw p ) < 
10 11 , show no sign of convergence at T > 10 4 c/w p (Keshet et al. 2009; Sironi 
et al. 2013). This is illustrated in Fig. 3 for a pair-plasma shock in 2D. 

For magnetized shocks, the situation is different, as we describe below 
(Sironi et al. 2013). At strong magnetizations, and for the quasi-perpendicular 
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field geometry most relevant for relativistic flows, particle acceleration is sup¬ 
pressed, and the shock quickly reaches a steady state. At low (but nonzero) 
quasi-perpendicular magnetization, the shock evolves at early times similarly 
to the case of unmagnetized shocks ( i.e. a = 0). Particle acceleration proceeds 
to higher and higher energies, and modes of longer and longer wavelength ap¬ 
pear. However, the maximum particle energy stops evolving once it reaches a 
threshold 7 sat oc cr -1 / 4 (Sironi et al. 2013), and at that point the overall shock 
structure approaches a steady state. (Sironi et al. 2013). 6 


4 PIC Simulations of Relativistic Shocks 

Only in the last few years, thanks to important advances in numerical algo¬ 
rithms and computer capabilities, plasma simulations have been able to tackle 
the problem of particle acceleration in relativistic shocks from first principles. 
In the following, we describe the major advances driven by large-scale PIC 
simulations in our understanding of particle acceleration in relativistic shocks. 
PIC codes can model astrophysical plasmas in the most fundamental way 
(Birdsall and Langdon 1991; Buneman 1993; Spitkovsky 2005), as a collection 
of charged macro-particles that are moved by the Lorentz force. The currents 
deposited by the macro-particles on the computational grid are then used to 
solve for the electromagnetic fields via Maxwell’s equations. The loop is closed 
self-consistently by extrapolating the fields to the macro-particle locations, 
where the Lorentz force is computed. 

Full PIC simulations can capture, from first principles, the acceleration 
physics of both electrons and ions. However, such simulations must resolve 
the electron plasma skin depth c/w pe , which is typically much smaller than 
astrophysical scales. Hence, most simulations can only cover limited time and 
length scales, and usually with low dimensionality (ID or 2D instead of 3D) 
and small ion-to-electron mass ratios (the ion skin depth c/ui vl is a factor of 
\Jrrii /m e larger than the electron skin depth c/w pe ). The results discussed 
below pertain to simulation durations of order ~ 10 3 — 10 4 w pe in electron- 
positron shocks and ~ 10 3 w p i in electron-ion shocks (but with reduced mass 
ratios), so a careful extrapolation is needed to bridge these microscopic scales 
with the macroscopic scales of astrophysical interest. Yet, as we review below, 
PIC simulations provide invaluable insight into the physics of particle injection 
and acceleration in astrophysical sources. 

The structure of relativistic shocks and the efficiency of particle acceler¬ 
ation depend on the conditions of the upstream flow, such as bulk velocity, 
magnetic field strength and field orientation. PIC simulations have shown that 
the shock physics and the efficiency of particle acceleration are insensitive to 

6 This conclusion regarding the saturation of the maximum particle Lorentz factor at 7 sa t 
has been tested in electron-positron shocks having a = 10 “ 4 — 10 —4 by Sironi et al. (2013), 
with the largest PIC study available to date. We caution that further nonlinear evolution, 
beyond the timespan covered by current PIC simulations, might be present in shocks with 
lower magnetization. 
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the shock Lorentz factor (modulo an overall shift in the energy scale), in the 
regime 1 of ultra-relativistic flows ( e.g. Sironi et al. 2013). Below, we 

only discuss results for shocks where the upstream Lorentz factor with respect 
to the downstream frame is q r > 5, neglecting the trans- and non-relativistic 
regimes that are outside the scope of this review. We discuss the physics of both 
electron-positron shocks and electron-ion shocks (up to realistic mass ratios), 
neglecting the case of electron-positron-ion shocks presented by e.g. Hoshino 
et al. (1992); Amato and Arons (2006); Stockem et al. (2012), which might 
be relevant for PWNe. As found by Sironi and Spitkovsky (2009a, 2011b); 
Sironi et al. (2013), for highly relativistic flows, the main parameter that con¬ 
trols the shock physics is the magnetization a. Below, we distinguish between 
shocks propagating into strongly magnetized media (er > 10~ 3 ) and weakly 
magnetized or unmagnetized shocks (a < 10 -3 ). 


4.1 Particle Acceleration in Strongly Magnetized Shocks 

For high magnetizations (a > 10~ 3 in electron-positron flows, or a > 3 x 10~ 5 
in electron-ion flows), the shock structure and acceleration properties depend 
critically on the inclination angle 9 between the upstream field and the shock 
direction of propagation (Sironi and Spitkovsky 2009a, 2011b). If the magnetic 
obliquity is larger than a critical angle 0 C rit, charged particles would need to 
move along the field faster than the speed of light in order to outrun the shock 
(“superluminal” configurations). In Fig. 4, we show how the critical angle 
$crit (as measured in the downstream frame) depends on the flow velocity 
and magnetization. In the limit of er <C 1 and q r 1, the critical obliquity 
approaches the value 0 C rit — 34°. 

Only “subluminal” shocks (9 < d cr ;t) are efficient particle accelerators 
(Sironi and Spitkovsky 2009a, 2011b; Sironi et al. 2013), in agreement with 
the analytical findings of Begelman and Kirk (1990). As illustrated in Fig. 5, 
a stream of shock-accelerated particles propagates ahead of the shock (panel 
(c)), and their counter-streaming with the incoming flow generates magnetic 
turbulence in the upstream region (panel (b)). In turn, such waves govern the 
acceleration process, by providing the turbulence required for the Fermi mech¬ 
anism. In the particular case of Fig. 5 — a relativistic shock with q r = 15, 
er = 0.1 and 9 = 15° propagating into an electron-ion plasma — the upstream 
turbulence is dominated by Bell-like modes (Reville et al. 2006; Lemoine and 
Pelletier 2010, 2011). The downstream particle spectrum in subluminal shocks 
shows a pronounced non-thermal tail of shock-accelerated particles with a 
power-law index 2 < s 7 < 3 (panel (d)). The tail contains ~ 5% of the parti¬ 
cles and ~ 20% of the flow energy at time 2250 w”; 1 ; both values appear to be 
time-converged, within the timespan covered by our simulations. 

In contrast, superluminal shocks {9 > 0 cr it) show negligible particle acceler¬ 
ation (Gallant et al. 1992; Hoshino 2008; Sironi and Spitkovsky 2009a, 2011b; 
Sironi et al. 2013). Here, due to the lack of significant self-generated turbu¬ 
lence, charged particles are forced to slide along the background field lines, 
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Fig. 4 Critical obliquity angle 0 t:T n (measured in the downstream frame) that separates 
subluminal and superluminal configurations (Sironi and Spitkovsky 2009a), as a function 
of the flow Lorentz factor 7 r and the magnetization tr, as indicated in the label. The filled 
black circle indicates our reference case with 7 r = 15 and a = 0.1. 


whose orientation prohibits repeated crossings of the shock. This inhibits the 
Fermi process, and in fact the particle distribution behind superluminal shocks 
is purely thermal. The same conclusion holds for both electron-positron and 
electron-ion flows. In electron-ion shocks, the incoming electrons are heated 
up to the ion energy, due to powerful electromagnetic waves emitted by the 
shock into the upstream medium, as a result of the synchrotron maser instabil¬ 
ity (studied analytically by Lyubarsky (2006), and with ID PIC simulations 
by e.g. Langdon et al. (1988); Gallant et al. (1992); Hoshino et al. (1992); 
Hoshino (2008)). Yet, such heating is not powerful enough to permit an effi¬ 
cient injection of electrons into the Fermi acceleration process at superluminal 
electron-ion shocks. 

If magnetized superluminal shocks are responsible for producing the radi¬ 
ating particles in astrophysical relativistic sources, the strong electron heating 
observed in electron-ion shocks implies that the putative power-law tail in the 
electron spectrum should start from energies higher than the ion bulk kinetic 
energy. For models of GRBs and AGN jets that require a power-law distri¬ 
bution extending down to lower energies, the presence of such shocks would 
suggest that electron-positron pairs may be a major component of the flow. 
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Fig. 5 Structure of an electron-ion subluminal shock with 7 r = 15, a = 0.1 and 6 = 15°, 
from Sironi and Spitkovsky (2011b). The simulation is performed in the downstream frame. 
The shock front is located at x ~ 725c/cj p i (vertical dotted red line in panel (a)), and it 
separates the upstream region (to its right) from the compressed downstream region (to its 
left). A stream of shock-accelerated ions propagates ahead of the shock (see the diffuse cloud 
in the momentum space x — p x i of panel (c) to the right of the shock, at x > 725 c/uj p i). 
Their interaction with the upstream flow (narrow beam to the right of the shock in panel 
(c)) generates magnetic turbulence ahead of the shock (see the transverse waves in panel (b), 
to the right of the shock). In turn, such waves govern the process of particle acceleration. 
In fact, the particle spectrum behind the shock (solid lines in panel (d); red for ions, blue 
for electrons) is not compatible with a simple thermal distribution (dashed lines), showing 
a clear non-thermal tail of high-energy particles, most notably for ions. 


4.2 Particle Acceleration in Weakly Magnetized and Unmagnetized Shocks 

Weakly magnetized shocks (a < 10~ 3 in electron-positron flows, a < 3 x 10 -5 
in electron-ion flows) are governed by electromagnetic plasma instabilities 
(see §3.1), that generate magnetic fields stronger than the background field. 
Such shocks do accelerate particles self-consistently, regardless of the mag¬ 
netic obliquity angle (Spitkovsky 2008a, b; Martins et al. 2009; Haugbpllc 2011; 
Sironi et al. 2013). The stream of shock-accelerated particles propagates ahead 
of the shock, triggering the Weibel instability. The instability generates fila¬ 
mentary magnetic structures in the upstream region, as shown in Fig. 6, which 
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Fig. 6 Shock structure from the 3D PIC simulation of a a = 10 electron-positron shock 
with 7 r . = 15, from Sironi et al. (2013). The simulation is performed in the downstream frame 
and the shock propagates along +x. We show the xy slice of the particle number density 
(normalized to the upstream density), and the xz and yz slices of the magnetic energy 
fraction eg. A stream of shock-accelerated particles propagates ahead of the shock, and their 
counter-streaming motion with respect to the incoming flow generates magnetic turbulence 
in the upstream via electromagnetic micro-instabilities. In turn, such waves provide the 
scattering required for particle acceleration. 


in turn scatter the particles back and forth across the shock, mediating Fermi 
acceleration. 

The accelerated particles in weakly magnetized shocks populate in the 
downstream region a power-law tail dN/dj oc 'y~ s x with a slope s 7 ~ 2.5, that 
contains ~ 3% of the particles and ~ 10% of the flow energy.' In electron-ion 
shocks, the acceleration process proceeds similarly for the two species, since 
the electrons enter the shock nearly in equipartition with the ions, as a result of 
strong pre-heating in the self-generated Weibel turbulence (Spitkovsky 2008a; 
Martins et al. 2009; Sironi et al. 2013). In both electron-positron and electron- 
ion shocks, the maximum energy of the accelerated particles scales in time 

as Tmacc oc f 1 / 2 (Sironi et al. 2013), as shown in Fig. 7. More precisely, the 

maximum particle Lorentz factor in the downstream frame scales as 

^jmax — 0.5 y r (cUp e f) / (4) 

7 max,i ~ ^ 0.25 7r Kit) 1 / 2 (5) 

in electron-positron and in electron-ion shocks, respectively (Sironi et al. 2013). 
This scaling is shallower than the so-called (and commonly assumed) Bohm 

7 These values are nearly independent of the flow composition and magnetization, in 
the regime of weakly magnetized shocks. Also, they are measured at time ~ 10 4 ujpe in 
electron-positron shocks and at ~ 10 3 o;~^ in electron-ion shocks, but they appear remark¬ 
ably constant over time, within the timespan covered by our simulations. 
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Fig. 7 Temporal evolution of the downstream particle spectrum, from the 2D simulation 
of a 7 r = 15 electron-ion (rrii/m e = 25) shock propagating into a flow with magnetization 
<j = 10 —5 , from Sironi et al. (2013). The evolution of the shock is followed from its birth 
(black curve) up to o; p i t = 2500 (red curve). In the top panel we show the ion spectrum and 
in the bottom panel the electron spectrum. The non-thermal tails approach at late times a 
power law with a slope s 1 = 3.0 for ions and s 1 = 2.5 for electrons (black dashed lines in 
the two panels). In the bottom panel, we overplot the ion spectrum at cj p i t = 2500 with a 
red dotted line, showing that ions and electrons are nearly in equipartition. Inset of the top 
panel: mean downstream ion (red) and electron (blue) energy, in units of the bulk energy 
of an upstream particle. The dashed blue line shows the electron energy at injection. Inset 
of the bottom panel: temporal evolution of the maximum Lorentz factor of ions (red) and 
electrons (blue), scaling as oc (cj p i£) 1//2 at late times (black dashed line). 


limit 7 max oc t , and it naturally results from the small-scale nature of the 
Weibel turbulence generated in the shock layer (see Fig. 6). 

The increase of the maximum particle energy over time proceeds up to 
a saturation Lorentz factor (once again, measured in the downstream frame) 
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that is constrained by the magnetization a of the upstream flow 

7 sat 




~ 2 j r a 1 ^ 4 


( 6 ) 

(7) 


in electron-positron and electron-ion shocks, respectively. The saturation of 
the maximum particle energy is shown in Fig. 8 for a shock with a = 10 -3 . 
Further energization is prevented by the fact that the self-generated turbulence 
is confined within a region of thickness Ls^at oc er^ 1 / 2 around the shock (Sironi 
et al. 2013). 



10° io 1 io 2 io 3 
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Fig. 8 Time evolution of the downstream particle spectrum from the 3D PIC simulation of 
a a = 10 -3 electron-positron shock with 7 ^ = 15, from Sironi et al. (2013). The evolution 
of the shock is followed from its birth (black curve) up to uj pe t = 3000 (red curve). We 
overplot the spectrum at u) pe t = 3000 from a 2D simulation with the same parameters 
(red dotted line), showing excellent agreement at high energies. The inset shows that the 
maximum particle Lorentz factor grows as 7 max oc t 1 ' 2 , before saturating at 7 sa t oc a~ 1 / 4 . 
The results are consistent between 2D (dotted) and 3D (solid). 


5 Astrophysical Implications 

5.1 Acceleration of Ultra-High Energy Cosmic Rays 

Relativistic shock waves have long been considered as prime candidates for the 
acceleration of cosmic rays to the highest energies observed, E ~ IO 20 eV. In¬ 
deed, a naive extrapolation of the acceleration time scale in the sub-relativistic 
regime (t acc ~ t s c&tt/Pu: with £ S catt the scattering timescale) suggests that rel¬ 
ativistic shocks ( i.e . f3 u ~ 1) accelerate particles on shorter time scales than 
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non-relativistic shocks (he. /3 U -C 1), at a given t s ca tt- For given radiative loss 
and escape time scales, this implies that relativistic shocks would be acceler¬ 
ating particles to much higher energies than non-relativistic shocks. However, 
the situation is more complex than it appears; in particular, in relativistic 
shock waves, t sca tt may be much larger than usually assumed. 

As mentioned repeatedly in the previous paragraphs, particle acceleration 
in the relativistic regime 7 U / 3 U 1 around a steady planar shock wave, oper¬ 

ates only if intense micro-turbulence has been excited in the shock precursor, 
as demonstrated analytically (Lemoine et al. 2006), by Monte Carlo simu¬ 
lations (Niemiec et al. 2006) and by PIC simulations (Sironi et al. 2013); 
consequences for the acceleration of particles to ultra-high energies have been 
discussed in several papers, e.g. by Pelletier et al. (2008), Lemoine and Wax- 
man (2009), Lemoine (2011), Eichler and Pohl (2011), Bykov et al. (2012); 
Sironi et al. (2013) or more recently by Reville and Bell (2014). 

Scattering in small-scale turbulence leads to a downstream residence time 
t sC att ~ r 2 /(Aasc), with tl the Larmor radius of the particle and A sb the 
coherence length scale of the turbulence. This implies that the (shock frame) 
acceleration timescale t acc grows quadratically with the energy, which fits well 
the result seen in PIC simulations that the maximum energy grows as the 
square root of time. In other words, as the particle energy grows, the accel¬ 
eration timescale departs more and more from the Bohm estimate, which is 
generally used to compute the maximum energy. Comparing for instance the 
acceleration timescale, which is at least equal to the above downstream res¬ 
idence time, with the dynamical timescale r/ 7 „ in the shock rest frame (r 
denoting the radius of the shock front in the upstream rest frame), one finds 
a maximum energy E max < e5Br , with 6B the strength of the 

turbulent field expressed in the shock frame; the above maximal energy has 
been written in the upstream (observer) frame. The factor in the brackets gen¬ 
erally takes very small values, because A sb ~ c/w p while r is a macroscopic 
length scale; this maximal energy is thus far below the so-called Hillas estimate 
eSBr, which corresponds to a Bohm estimate for f scatt . 

Another way to phrase the problem is as follows (see Lemoine and Wax- 
man 2009 for a discussion): assume that the acceleration timescale is writ¬ 
ten f acc = Ai'l/c , and derive the maximum energy by comparing f acc with 
fdyn = r/( 7 /3c) as above for a jet moving at velocity /3 towards the observer. 
Then one finds that acceleration of particles of charge Z to 10 2 O l? 2 o eV re¬ 
quires that the isotropic equivalent magnetic luminosity of the object exceeds: 
Lb > 10 45 Z~ 2 E% 0 A 2 "/ 2 erg/s, a very large number indeed, all the more so 
if A 1. For acceleration at ultra-relativistic shock waves, A is much larger 
than unity (while the Bohm estimate corresponds to A ~ 1), with typical val¬ 
ues A ~ E/ ( 7 u m p c 2 ). In summary, particle acceleration at ultra-relativistic 
shock waves does not appear fast enough to produce particles of ultra-high 
energies. In particular, when the above arguments are applied to the case of 
the external shock of a GRB, the maximal energy is found to be of the order 
of 10 16 eV (Plotnikov et al. 2013; Sironi et al. 2013; Reville and Bell 2014). 
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It is important however to note three caveats in the above arguments. One 
is that as 7„/3 u —> 1 , i.e. for mildly relativistic shock waves, the nature of the 
turbulence remains unknown and one cannot exclude that scattering would 
be closer to a Bohm estimate. Two facts support such a speculation: (1) the 
precursor increases in size as ") u diminishes, which suggests that MHD-scale 
instabilities could arise and excite large scale turbulence; and (2), the obliquity 
becomes less of a problem for mildly relativistic shock waves, suggesting that 
large scale turbulence could possibly lead to acceleration in this regime. A sec¬ 
ond caveat is the fact that PWNe are very efficient particle accelerators, even 
though one would expect the opposite in the absence of reconnection or other 
dissipative processes, due to the large magnetization of the flow (Section 5.3). 
More precisely, synchrotron photons are observed with energies as high as 
100 MeV, which means that pairs are accelerated up to the radiation-reaction 
limit, i.e. with an acceleration time scale close to the theoretical Bohm scaling. 
Such empirical evidence suggests that ions could also be accelerated to very 
high energies, if ions are indeed injected along with pairs in the wind. In the 
Crab, such a maximal energy would be limited by the confinement in the neb¬ 
ular turbulence to values of the order of 10 1 ' eV (for Z = 1); more powerful 
nebulae, associated with young pulsars born with a few millisecond periods, 
could however confine (and potentially accelerate) protons up to the highest 
energies (Lemoine et al. 2014). Finally, as the nonlinear evolution of weakly 
magnetized or parallel shocks over long timescales is not yet understood, some 
of the above estimates, pertaining e.g. to the diffusive properties and extent 
of the magnetic field, may be altered on macroscopic times. 


5.2 Radiative Signatures of Relativistic Blast Waves 

In line with the previous discussion, one can compute the maximal energy for 
electrons and derive the maximal synchrotron photon energy. Using an acceler¬ 
ation time scale t acc ~ r 2 L /(Xsbc) and comparing to synchrotron losses in the 
self-generated micro-turbulence, characterized by its magnetization es, one de¬ 
rives a maximum synchrotron photon energy of the order of a few GeV in the 
early phase of GRB afterglows, i.e. during the first hundreds of seconds (Kirk 
and Reville 2010; Plotnikov et al. 2013; Lemoine 2013; Wang et al. 2013; Sironi 
et al. 2013). Let us stress that in the latter study, this estimate has been derived 
from PIC simulations with a self-consistent measurement of the acceleration 
time scale in the self-consistent magnetic field. The synchrotron radiation of 
electrons accelerated at the external ultra-relativistic shock of GRBs can thus 
produce the bulk of the long-lasting > 100 MeV emission detected by the 
Fermi satellite (e.g. Kumar and Barniol Duran 2009; Ackermann 2010; De 
Pasquale 2010; Ghisellini et al. 2010). The photons that have been observed 
with energies in excess of > 10 GeV probably result from inverse Compton 
interactions (Wang et al. 2013). Interestingly, the recent GRB130427A has 
revealed a long-lasting emission with a possible break in the spectrum at an 
energy of a GeV, characteristic of a turn-over between the synchrotron and the 
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synchrotron self-Compton components (Tam et al. 2013), in good qualitative 
agreement with the above arguments. 

Other potential radiative signatures of the shock microphysics come from 
the small-scale nature of the turbulence and its long-term evolution in the 
blast. As discussed in Section 3.2, one notably expects this turbulence to relax 
through collisionless damping on hundreds of c/w p (Chang et al. 2008; Keshet 
et al. 2009; Lemoine 2015) while the electrons typically cool on much longer 
length scales. In CRB external blast waves, the shocked region is typically 7—9 
orders of magnitude larger than c/w p in size, which leaves room for a substan¬ 
tial evolution of e^, even if it decreases as a mild power-law in distance from 
the shock, as suggested by the above studies. Since the electron cooling length 
depends on the inverse of the electron Lorentz factor, particles of different 
initial Lorentz factors emit their energy in regions of different magnetic field 
strength, leading to a non-standard synchrotron spectrum (Rossi and Rees 
2003; Derishev 2007; Lemoine 2013), which could in principle be used as a 
tomograph of the evolution of the micro-turbulence downstream of the shock. 
Interestingly, in this picture the decay index of the turbulence is related to 
the long-wavelength content of the power spectrum of magnetic fluctuations 
at the shock front, which is unknown so far, as it is known to be modified by 
the acceleration of higher energy particles (Keshet et al. 2009). Finally, it is 
interesting to note that the recent broad-band analysis of GRB afterglows seen 
from the radio up to GeV energies has indeed revealed spectral signatures of 
a decaying magnetic field (Lemoine et al. 2013), with a decay law scaling with 
distance from the shock roughly as At -0,5 (Ax being the proper distance to 
the shock in the downstream frame). 

As discussed in Section 3.2, there are alternative possibilities however; it 
has been suggested for instance that the turbulence could evolve in a self- 
similar way as a function of distance to the shock, maintaining a uniform es 
thanks to an inverse cascade process (Katz et al. 2007). It is also possible that 
external sources seed the blast with a large scale long-lived turbulence, e.g. 
through a Rayleigh-Taylor instability at the contact discontinuity (Levinson 
2009) or through small scale dynamos following the interaction of the shock 
front with external inhomogeneities (Sironi and Goodman 2007; Couch et al. 
2008). Hopefully, future high accuracy observational data will provide diag¬ 
nostics which can be confronted with numerical simulations. 

The possibility that the small scale nature of the turbulence gives rise to 
diffusive (or jitter) synchrotron radiation rather than conventional synchrotron 
radiation has also attracted attention (e.g. Medvedev 2000, 2006; Fleishman 
2006; Mao and Wang 2011; Medvedev et al. 2011; Kelner et al. 2013). In par¬ 
ticular, jitter radiation has been proposed as a solution for the fact that GRBs 
prompt spectra below the peak frequency are not always compatible with the 
predictions of synchrotron emission (the so-called “line of death” puzzle, see 
Preece (1998)). In the jitter regime, particles are deflected by less than l/q 
(7 is the electron Lorentz factor) as they cross a wavelength A sb, implying 
that coherence of the emission is maintained over several coherence cells of 
the turbulence. This regime thus takes place whenever the wiggler parame- 
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Fig. 9 Ab initio photon spectrum (thick solid lines) from the 2D PIC simulation of an 
unmagnetized (i.e. a = 0) pair shock. Red lines are for head-on emission (h = x, along 
the shock direction of propagation), blue lines for edge-on emission (h = y. along the shock 
front). The slope at low frequencies is 2/3 (black long-dashed lines), proving that the spec¬ 
tra are consistent with synchrotron radiation from a 2D particle distribution (in 3D, the 
predicted slope of 1/3 is obtained). By separating the relative contribution of downstream 
( x < .r s h ; thin solid lines) and upstream (x > .x H [,: dotted lines) particles, one sees that up¬ 
stream particles contribute significantly to the total emission (thick solid lines), especially 
at high frequencies. Frequencies are in units of the plasma frequency tu p . 


ter a = eSBXgg/mcr <C 1, while the standard synchrotron approximation 
becomes valid in the opposite limit. However, it is easy to verify that in the 
vicinity of the shock a ~ 7 | sh , with 7 | sh the average Lorentz factor of the 
supra-thermal electrons in the shock rest frame, suggesting that jitter signa¬ 
tures must be weak. 

The absence of jitter radiation in relativistic shocks has been demonstrated 
from first principles by computing the radiation from particles in PIC simula¬ 
tions (Sironi and Spitkovsky 2009b) , which produce spectra entirely consistent 
with synchrotron radiation in the fields generated by the Weibel instability 
(Fig. 9). The so-called “jitter” regime is recovered only by artificially reducing 
the strength of the fields, such that the parameter a becomes much smaller 
than unity. So, if the GRB prompt emission results from relativistic unmagne- 
tized shocks, it seems that resorting to the jitter regime is not a viable solution 
for the “line of death” puzzle. At frequencies above the peak, the synthetic 
spectra from PIC simulations show, somewhat unexpectedly, that the contri¬ 
bution of the upstream medium to the total emission is not negligible (Fig. 9), 
yet it is omitted in most models. This causes the radiation spectrum to be 
flatter than the corresponding downstream spectrum, thus partly masking the 
contribution of downstream thermal particles. 
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5.3 Radiative Signatures of Pulsar Wind Nebulae 

The spectrum of PWNe consists of two components, where the low energy 
component, most likely dominated by synchrotron, shows a cutoff at a few tens 
of MeV. The fact that synchrotron emission reaches these energies, despite the 
rapid synchrotron cooling, implies that particle acceleration in the nebula is 
an extremely fast process (de Jager and Harding 1992), which challenges our 
understanding of particle acceleration in relativistic shocks. 

Around the equatorial plane of obliquely-rotating pulsars, the wind con¬ 
sists of toroidal stripes of opposite magnetic polarity, separated by current 
sheets of hot plasma. It is still a subject of active research whether the al¬ 
ternating stripes will dissipate their energy into particle heat ahead of the 
termination shock, or whether the wind remains dominated by Poynting flux 
till the termination shock (Lyubarsky and Kirk 2001; Kirk and Skjaeraasen 
2003; Sironi and Spitkovsky 2011a). If the stripes are dissipated far ahead 
of the termination shock, the upstream flow is weakly magnetized and the 
pulsar wind reaches a terminal Lorentz factor (in the frame of the nebula) 
7 r ~ L sd /m e c 2 N ~ 3.7xl0 4 L sdj38 .5A r ^o 1 , where L sd = 3x10 3S L sd , 3 $, 5 ergss -1 
is the spin-down luminosity of the Crab (the Crab Nebula is the prototype of 
PWNe), and N = 10 4o A4o s~ 4 is the particle flux entering the nebula, includ¬ 
ing the radio-emitting electrons (Bucciantini et al. 2011). 

For electron-positron flows, as appropriate for pulsar winds, the maxi¬ 
mum particle Lorentz factor in the downstream frame increases with time 
as "fmax ~ 0.5 7 r (wpi) 1 / 2 (see Section 4). The plasma frequency u> p can be 
computed from the number density ahead of the termination shock, which 
is uts = iV/(47ri?Tgc), assuming an isotropic particle flux. Here, Rts = 
3 x 10 17 i?TS,i7.5 cm is the termination shock radius. Balancing the acceleration 
rate with the synchrotron cooling rate in the self-generated Weibel fields, the 
maximum electron Lorentz factor is 



( 8 ) 


7 sync,e — 8.5 X 10 L 


A stronger constraint comes from the requirement that the diffusion length 
of the highest energy electrons be smaller than the termination shock radius 
( i.e . a confinement constraint). Alternatively, the acceleration time should be 
shorter than i?Ts/c, which yields the critical limit 



(9) 


which is generally more constraining than the cooling-limited Lorentz factor 
7sj/nc,e- The corresponding synchrotron photons will have energies 


hvconf ,e — 0*17 L 



( 10 ) 


which are apparently too small to explain the X-ray spectrum of the Crab, 
extending to energies beyond a few tens of MeV. We conclude that Fermi 
acceleration at the termination shock of PWNe is not a likely candidate for 
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producing X-ray photons via the synchrotron process, and valid alternatives 
should be investigated. 

One possibility - magnetic dissipation of the striped pulsar wind in and 
around the shock front itself has been extensively studied, with the conclu¬ 
sion that particle acceleration along extended X-lines formed by tearing of the 
current sheets may contribute to the flat particle distribution (with spectral 
index s 7 ~ 1.5) required to explain the far infrared and radio spectra of PWNe 
( e.g. , Lyubarsky 2003; Sironi and Spitkovsky 2011a, 2012). Indeed, hard par¬ 
ticle spectra are found to be a generic by-product of magnetic reconnection 
in the relativistic regime appropriate for pulsar winds (Sironi and Spitkovsky 
2014; Sironi et al. 2015, see also Kagan et al. (2015) in the present volume). 
However, further acceleration to gamma-ray emitting energies by the Fermi 
process cannot occur in the shock that terminates the pulsar wind, if particle 
scattering depends only on the self-generated turbulence. 

Yet, the steady-state hard X-ray and gamma-ray spectra of PWNe do 
look like the consequences of Fermi acceleration - particle distributions with 
s 7 ~ 2.4 are implied by the observations. In this regard, we argue that the 
wind termination shock might form in a macroscopically turbulent medium, 
with the outer scale of the turbulence driven by the large-scale shear flows in 
the nebula (Komissarov and Lyubarsky 2004; Del Zanna et al. 2004; Camus 
et al. 2009). If the large-scale motions drive a turbulent cascade to shorter 
wavelengths, back-scattering of the particles in this downstream turbulence, 
along with upstream reflection by the transverse magnetic field of the wind, 
might sustain Fermi acceleration to higher energies. 

Another “external” influence of reconnection on the shock structure, that 
might lead to particle acceleration to higher energies, may be connected to 
the accelerator behind the recently discovered gamma-ray flares in the Crab 
Nebula (Abdo 2011). Runaway acceleration of electrons and positrons at re¬ 
connection X-lines, a linear accelerator, may inject energetic beams into the 
shock, with the mean energy per particle approaching the whole open field 
line voltage, > 10 16 V in the Crab (Arons 2012), as required to explain the 
Crab GeV flares. This high-energy population can drive cyclotron turbulence 
when gyrating in the shock-compressed fields, and resonant absorption of the 
cyclotron harmonics can accelerate the electron-positron pairs in a broad spec¬ 
trum, with maximum energy again comparable to the whole open field line 
voltage (Hoshino et al. 1992; Amato and Arons 2006). 


6 Conclusions 

There has been significant progress in our understanding of relativistic shocks 
in recent years, thanks to both analytical work and numerical simulations. The 
highly nonlinear problem of particle acceleration and magnetic held generation 
in shocks — with the accelerated particles generating the turbulence that in 
turn mediates their acceleration — is being tackled from first principles, as¬ 
sessing the parameter regime where particle acceleration in relativistic shocks 
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is efficient. In this chapter, we have described the basic analytical formalism 
of test particle acceleration in relativistic shocks, leading to the “universal” 
energy slope s 7 ~ 2.2 in the ultra-relativistic limit; we have unveiled the most 
relevant plasma instabilities that mediate injection and acceleration in rela¬ 
tivistic shocks; and we have summarized recent results of large-scale PIC simu¬ 
lations concerning the efficiency and rate of particle acceleration in relativistic 
shocks, and the long-term evolution of the self-generated magnetic turbulence. 
Our novel understanding of particle acceleration and magnetic field generation 
in relativistic shocks has profound implications for the modeling of relativistic 
astrophysical sources, most importantly PWNe, GRBs, and AGN jets. 
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